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Abstract 



J> ' The multiplicity distributions produced by the variation of time-dependent gravitational fields 

in a conformally flat background geometry belong to the same class of infinitely divisible distribu- 
^T) . tions found, for fixed centre of mass energies and symmetric (pseudo)rapidity intervals, in charged 

multiplicities produced in pp, pp and in heavy ion collisions. Apparently unrelated multiplicity 
distributions are classified in terms of the (positive) discrete representations of the SU(1, 1) group. 
The gravitational analogy suggest a global high-energy asymptote for the distributions measured 
in pp and pp collisions. Second-order cross correlations between positively and negatively charged 



distributions represent a relevant diagnostic for a closer scrutiny of the multiparticle final state. 
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A. Sakharov [1] was presumably the first to raise the question of the quantum mechanical origin of 
density perturbations in the early Universe suggesting that the complicated patterns observed in the 
galaxy distributions could have a plausible origin in the zero-point fluctuations of matter and radiation 
fields in curved backgrounds. More recently the latter perspective gained a firmer support from the 
analyses of the Cosmic Microwave Background (CMB) anisotropies and polarization. It can therefore 
be speculated that the initial conditions of the CMB anisotropies were actually set in a regime where 
the expansion rate was H ~ O(10~ 6 Mp) ~ O(10 37 Hz). Today the background geometry is, to a good 
approximation, conformally flat at least according to the WMAP data in their various releases [2] (see 
also, for instance, [3]) and in the light of the so-called ACDM lore where A stands for the cosmological 
constant and CDM for the cold dark matter component. 

Gravitational and strong interactions are intrinsically different in the physical regimes where direct 
laboratory tests are possible. Even acknowledging, for large space-time curvatures, the corresponding 
largeness of the expansion rate, multiple production in strong gravitational fields is naively expected to 
be drastically different from the laws governing the dynamics of the multiparticle final state in hadronic 
processes such as, for instance, proton-proton (i.e. pp), proton-antiproton (i.e. pp) or even heavy ions 
collisions. The purpose of the present analysis is to challenge the latter statement by analyzing the 
multiple production of electrically charged species (i.e. positive and negative) in conformally flat 
geometries and in a correct dynamical framework. Consider therefore one of the simplest models for 
the production of charged species, i.e. a minimally coupled complex scalar field in a conformally flat 
geometry; the relevant four-dimensional action can be written 



S = / d 4 x \/—g 



g^da<l>*d v <t> - m 2 </>> 



(1) 



where g^ u = a 2 (T)r]^ v is the (conformally flat) background metric of Friedmann- Robertson- Walker 
type expressed in terms of the conformal time coordinate r and in terms of the scale factor a{r)\ 
rj^ u denotes the Minkowski metric with signature (+, — , — , — ). By introducing the canonical normal 
modes &(x, r) = a(r)(f>(x,T) and <&*(x, r) = a(r)(p*(x,T) the canonical Hamiltonian can be written as 



H(t) = / d 3 x 



n*n + h($ n + $*n*) + d^d 1 ® + m V$*$ 



(2) 



where % = d T (In a) denotes the derivative of the natural logarithm of the scale factor with respect 
to the conformal time coordinate and n = d T <&* — "H<1>* is the canonical momentum. By promoting 
the classical fields to quantum mechanical operators obeying (equal time) commutation relations the 
Hamiltonian (2) can be written as the sum of a free part and of an interacting part H(t) = Hq(t) + 
Hj(t) with Hq(t) and -ffi(r) given, respectively, by 

H (t) = 2Jd 3 p w (r) K (p), H^t) = 2 J d 3 p [X*(r)K_(p) + X(r)K + (p)}, (3) 

where w(r) = \Jp 2 + m 2 a 2 , pis the comoving three-momentum and A(r) = iTi/2; the operators K±(p) 



and Kq(p) are denned as: 



k +(p) = Opbip, K_(p) = apb-p, K (p) = - 



ajfOp + b-pblp 



(4) 



Since [ap, at] = 5^ 3 >(p — k), [b$, bt] = 5^\p — k) and [ag, br] = 0, the operators (4) satisfy the 
commutation relations 

[K_(p),K + (q)] = 2K (p)5^(p- q), [K (p),K±(q)] = ±K ± (p)5^(p- q), (5) 

showing that K±(p) and Kq{p) are nothing but the generators the SU(1, 1) group obeying the com- 
mutation relations of the corresponding Lie algebra. Owing to the group structure (5) and to the 
specific form of the Hamiltonian of Eq. (3), the multiparticle state for r — > +oo can be obtained by 
applying to the initial state \^i(p)) the product of two unitary operators lZ((p p ) and Ti(z p ): 

\* f (p)) = K(ip p ) S(z p )|^(p)), |*/) = II !*/(#>> ( 6 ) 

p 

where \^>f(p)) denotes the final state and where the unitary operators are given by: 

K(ip p ) = exp \-2itpp K (p)], T,(z p ) = exp [z* K-(p) - z p K + (p)], (7) 

with z p = r p e l ® p and a v = (2ip p — $„). Denoting with a an overdot a derivation with respect to 
r, the time evolution of the variables r p (r), <£> p (t) and a p (r) is given by 2 r p = —7-Lcosa p , ip p = 
uj + % tanh r p sin a p , and a p = 2oj + 2% sin a p / tanh 2r p . The explicit solution of the latter equations 
depends upon the evolution of H which is determined, ultimately, by the evolution of the space-time 
curvature. The purpose here will not be to compute the average multiplicity in a given model but 
rather to analyze the correlation properties of the multiplicity distributions. 

The initial quantum state for each mode of the field (i.e. |vPj(p))) can be classified in terms of the 
Fock basis |n+ n_) which is an appropriate basis for the irreducible representations of SU(1, 1) once the 
group generators are represented as in Eq. (4) (see also [8] for the SU{2) case). An equivalent basis for 
the irreducible representations of SU{1, 1) is provided by the vectors \qn<±) where q = n + — n_ is the 
total charge and n c h = n + + n^ is the total number of charged species. The vectors \qn c ^) are the stan- 
dard basis of the irreducible representations T + of SU(1, 1) where k is the principal quantum number 
and m is the magnetic quantum number, i.e. the eigenvalue of Kq. The negative series T~ k is sym- 
metric under the exchange n+ — > n_ while the principal (continuous) series will not play a specific role 
in the present considerations. In terms of k and m we have that the total charge and the total number 
of particles are given, respectively, by q = 2k — 1 and by n c h = 2m — 1. Since 2k = N is a positive in- 
teger, the charge and the total number of charged species will be, respectively, q = M — 1 = 0, 1, 2, ... 



2 It should be noted that the factorization of the time-evolution operator in terms of lZ((p p ) and of £(z P ) is non-trivial 
since Ho and Hi do not commute. An analog problem arises in two-photon optics [6, 7]. 



and n c h = q, q + 1, q + 2... and so on. According to the Backer-Campbell-Hausdorff (BCH) de- 
composition, the operator of Eq. (7) can be factorized as the product of the exponentials of the 
group generators [6, 9], i.e. £(z p ) = A+(z p )Ao(z p )A-(z p ) where Ao(z p ) = exp [— 21n (coshr p )Ko(p)], 
A-(zp) = exp [e~ l ® p tanh r p K-(p)], and A+(z p ) = exp [— e l ® p tanh r p K + (p)] . Using the BCH decom- 
position, as well as the explicit form of \^f(p)) given in Eq. (6), the density operator can be deduced. 
Defining the following triplet of functions 

C(m £ i) ~ nm + l)T(m-l + j + l) 

M(n±,£,j,£',j') = Jc(n + ,£,J)C(n-,ej)C(n + ,e',j')C(n-,£',f), (9) 

Fir a -n+-£ £'■ i i') - e -'<M(J-l)-G <'-t ')] (-tannr p j 

j-{r p , a p ,n±, i, i ,j,j ) -e ^ ^ ^ n++n _ _ e _ e , +1 , {*■") 

the density operator reads 

PI(P)= E E G(n ± ;£,£ , ;j,j')\n + -£ + j,n^-£ + j)(j'-£' + n.,j-£ + n + \, (11) 
j,j'=o e,e=o 

a(n±; £,£'; j,j')=F(r p ,a p ;n±;£,£';j,j / )M(n±,£,j,£ / ,j'), (12) 

where £ max = Min(n + , n_). In terms of the density matrix of Eq. (11) the expectation value of a 
generic operator O can be computed as (O) = Tr[p O]. Consider first, in Eq. (11), the case q = n^ = 0. 
Since, in this case, ^ max = 0, Eq. (11) implies: 



PM = =— 1 E (-ir +n e-^ m - n H^i) \mm){nn\, (13) 



where (N+ + N_) = 2sinh r p = 2n and n denotes the average multiplicity of pairs. The density 
matrix of Eq. (13) is idempotent (i.e. p\ = pf) but, nonetheless the diagonal elements of pf(p) follow 
a Bose-Einstein (BE) multiplicity distribution: 

— n °° 

P^\n) = (nn\p f (p)\nn) = W , £ p( BE \n) = 1. (14) 

^ ' n=0 

The integration of Eq. (13) over a p between and 2ir defines the reduced density operator in the 
random phase approximation: the mixed state obtained with this procedure will have a thermal 
density matrix. This simple observation shows that a BE multiplicity distribution is not sufficient 
to infer local thermal equilibrium: as an example recall that chaotic light distributed as BE can 
be generated by sources in which atoms are kept at an excitation level higher than that in thermal 
equilibrium [10]. The multiplicity distributions of Eqs. (13) and (14) are also typical of the two-mode 
squeezed vacuum states [6, 11, 12]. The two-mode squeezed vacuum state of photons can be used to 
describe multigraviton [13], multiphonon [14], and multiphoton [15] states in the quantum treatment 
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Figure 1: The distributions of Eqs. (15) and (17) are compared for the same range of parameters. In 
both plots the scale is linear on the horizontal axis, while, on the vertical axis the common logarithm 
is reported. In the legends of both plots (n) coincides with the n appearing in Eqs. (15) and (17). 
The naughts on the two (lower) thin lines are not experimental points but just a graphic means to 
distinguish different curves reminding visually that we are dealing here with discrete distributions. 



of cosmological inhomogeneities (see also [16] for an introductory discussion). Consider next the case 
q = n c h = N — 1, i.e. the situation in which the initial state has a definite charge. Equation (11) 
implies that the diagonal elements of the density operator can be written as 



Pn{n,M) 



T(n + N) 



n 

AT 



r(A0r(n + 1) 



AT 



+ 1 



HAP 



■A/-1 
2M ' 



(15) 



which reduces to the second expression of Eq. (14) when M — Y 1. The multiplicity distribution of Eq. 
(15) falls into the class of negative binomial multiplicity distributions. Indeed, by expanding Eq. (15) 
in powers of e = (TV - l)/(2/V) < 1: 



P n {n,M)=P^ B \n,M) 



1 + (n — n)( — 



N 



Pn(n,k^ B ) 



n 



n+N 

n 



^)m-xX 

n 



k 



NB 



r(n + km) 

r(&NB)r(n + 1) \n + &NB/ V^ + ^NB 



&NB 



(16) 
(17) 



where P n {n, &nb) denotes the negative binomial (NB) probability distribution with parameters n and 
&nb- Equation (15) is in fact a distorted negative binomial distribution (DNB): according to Eq. (16) 
P n (n,J\f) is larger than P n (n,J\P) for n < n and it is smaller than P n (n,J\f) for n > n. In Fig. 1 
the DNB distribution of Eq. (15) is compared, for the same range of parameters, with the standard 
form of the NB distribution (see Eq. (17)). The full lines denotes the standard negative binomial 



result while the dashed lines refer to the DNB of Eq. (15) which undershoots (in comparison with 
the standard NB case) when the multiplicity exceeds n and overshoots in the opposite case, i.e. when 
n <n. In the plot at the right of Fig. 1 the notation is the same but the distributions are illustrated 
in Koba-Nielsen-Olesen (KNO for short) variables [17]. While the KNO scaling is violated since /cnb 
changes with the centre of mass energy of the collision (see below), there is still the useful habit 
to present the results for the multiplicity distributions in terms of KNO variables. To complete the 
discussion of the limits of Eq. (11) it is interesting to mention, the case when n c h ^ but q = 0; the 
density matrix of Eq. (11) leads then to a more complicated multiplicity distribution whose explicit 
form can be written as 

1 I „\t 12 



P m (n,x) 



r rV + i)r"(~ + 1) [g r(m _ n + t + ^ {( ' + 1)r3(n 



(x + l)m+n-l 

sinh r p = (2n — n)/[2(n + 1)]. In the case n = we have that 



+ D 



(18) 



Pm with X 



0: the finite 
n. The 



where x 

sum appearing in Eq. (18) gives, after squaring, 1/T 2 (m + 1) and P m (x) 

multiplicity distributions described by Eq. (18) have several interesting features like oscillations in m 

for fixed values of ft and n. These properties will not be specifically discussed here. Neglecting some 

phases, the connection between the multiplicity distributions and the SU(1, 1) group structure can be 

neatly expressed by computing, in explicit terms, the Wigner matrix element of the positive discrete 

scries: 



^+k 



{km'\T 1 {z p )\km) 
F[l — m — k, k - 



y/T{k + m') T(m' - k + 1) 



i l ®v tanh r„ 



y/T(k + m) T(m-k + l) (m' - m)\ (coshr p ) 2m 



m; m 



m + 1; — sinh r 



Ph 



m > m, 



where F[a, (3;y;x] is the hyper geometric function. From Eq. (19) we have, immediately, that 

2 



T l/2 
1/2 n+1/2 

e/2 n+i/2 



1 



n 



n + 1 V" + 1 
T(n + 1)T(£) \n + 



p(BE) ( . 



nj 



n 



sinh r. 



p- 



+ n 



pr\nj), 



sinh r r , 



n 



(19) 

(20) 
(21) 



The multiplicity distribution P n (n, /cnb) has various notable limits. Indeed, by computing the corre- 
sponding probability generating function (pgf ) , 



00 fc NB 

V(z,n,k NB ) = ^z k P k (n,k NB ) = { _ {i _ z) + feNB]fcNB , 



(22) 



it is immediate to see that, in the limit /cnb — >■ 1 we recover the pgf of the Bose-Einstein distribution 
while, in the limit &nb - > o° (and n fixed) V(z,n, /cnb) — > exp [(z — l)n], i.e. the pgf of the Poisson 
distribution. Finally, both the Gamma as well as the logarithmic distributions can be obtained as 
specific limits of the NB distribution. The distributions arising as matrix elements of SU(1, 1) from 



the recipe of Eqs. (20) and (21) are all infinitely divisible 3 . While a more thorough investigation of 
the latter statement is beyond the scopes of the present discussion, it is suggestive to note that neither 
the binomial distribution nor the uniform distribution (which are not infinitely divisible) arise in the 
proposed correspondence between SU(1, 1) matrix elements and the discrete multiplicity distributions. 

Multiplicity distributions in general (and the gravitational multiplicity distributions in particular) 
can be classified by using the degree of second order coherence. Borrowing the terminology of quantum 
optics, multiplicity distributions can be classified according to the their normalized two-point function 

[10, 11] 

(2) , _ (alaJapap,) ~ (N + ) 2 (2) {bj/bl/b_ f b_ p ) - (N-) 2 

.'/+ • ~^—, • !J- i " 7^-p , {^) 



<44 


Opttp", J 


(N + ) 2 


<4^ 


(N + ) 2 


> 
-{N + ){N_) 



(2) _ 1 = , v - v , ^,, __j ■__ (24) 

(N-)(N + ) 

(2) (2) (2) 

If q = n c h = 0, Eqs. (23) and (24) imply g K + ' = g_ = 2 while g± = 2 + 1/n, where, as in Eq. 

(14), n = sinh r p denotes the averaged multiplicity of pairs. If q = but n c h 7= 0, g + = g_ 7^ g± . 

(2) (2) (2) 
Finally, if q = n c h / we shall have that g\ 7^ g_ 7^ g± . To infer the correlation properties of 

the multiparticle final state the multiplicity distribution of the positive and negative species should 

(2 
be separately assessed and, eventually, cross-correlated to measure g± . 

The multiplicity distributions measured in hadronic reactions or even heavy ion collisions more 
often than not count all the charged species in the final state (not only the positively or negatively 
charged species). In analogy with the gravitational case, second-order (cross) correlations between 
positively and negatively charged distributions may represent a valid tool for the scrutiny of the 
multiparticle final state. Absent the latter measurements, and only for purposes of comparison, it is 
therefore simpler to consider a single mode of the field and to recall that 

m D 2 -IN) 
< g( 2 > - 1 = ^—L < 1, 25 

(iV}2 " V ' 

where D 2 = (N 2 ) — (N) 2 is the dispersion. In Eq. (25) g^ — 1 is defined as in Eq. (23). The 
expression in terms of the variance D 2 follows by recalling that (a* a' aa) = (N 2 ) — (N) where 
TV = a' a. If g( 2 > = 1, D 2 = (N) as it happens for a standard coherent state characterized by a Poisson 
multiplicity distribution. According to the quantum optical terminology, if g^ 2 ' > 1 the light is said to 



3 A probability distribution is said infinitely divisible if, for any given non negative integer k, it is possible to find k 
independent identically distributed random variables whose probability distributions sum up to the original distribution. 
The generating function of the sum of independent identically distributed random variables is given by the product of the 
generating functions of each distribution of the sum. On the basis of the latter theorem, a distribution with probability 
generating function V(z) is infinitely divisible provided, for any integer k, there exist k independent identically distributed 
random variables with generating function Qk(z) such that V(z) = [Qk(z)] k . 



be bunched (with super-Poissonian field statistics) while if g^ 2 ' < 1 the light is said to be antibunched 
(with sub-Poissonian field statistics) [10, 11]. In the case of a (single mode) Fock state \n) it can be 
readily shown from the definition (25) that g^ 2 > = 1 — 1/n. In more general terms the degree of second 
order coherence of an arbitrary single-mode excitation must satisfy g( 2 ' > 1 — \/{N) where the equality 
is saturated in the case of a Fock state. In short the gravitational multiplicity distributions derived in 
Eqs. (14), (15)— (17) and (18) are characterized by pronounced tails exhibiting a degree of correlation 
which can be quantitatively assessed from the degree of second-order coherence g^ 2 > ranging between 
1 and 2. 

The gravitational multiplicity distributions will now be compared and contrasted with the analog 
observables typical of hadronic collisions (see [18, 19, 20, 21] for a collection of relevant review articles 
on multiple production in strong interactions). Hadronic collisions at high energies lead to charged 
multiplicity distributions whose shapes are well fitted by a single negative binomial distribution in 
fixed intervals of central (pseudo) rapidity r/ [23]. There are of course physical differences between the 
charged multiplicity distributions arising from purely leptonic initial states (such as electron-positron 
collisions) and purely hadronic initial states such as {pp and pp). The physical analysis of different 
reactions in a unified dynamical framework (such as, for instance, the clan model [22]) would be 
interesting per se and it is anyway beyond the aims of the present investigation. In what follows 
the attention will be limited, for practical reasons, to the recent pp results at the LHC [23, 24, 25]. 
Denoting with the prime a derivation with respect to z of the probability generating function of Eq. 
(22), the mean and the dispersion of the distribution are given, respectively, by n = V(l) and by 
D 2 = V"(l) + V'{1) - [V'{1)} 2 and, consequently, 

5 (2) -l = 7^-A = I J -' C 2 = g^ + L (26) 

(n) z (n) fe NB n 

where C q = (n q )/(n) q are the normalized moments of the distribution and, in particular, C 2 = 
(n 2 }/(n) 2 ; the experimental collaborations instead of reporting the values of &inb prefer to report 
sometimes the value of C 2 (see, e.g. second and third paper of [23]). It is a property of the negative 
binomial distribution that all the higher-order moments (i.e. C q with q > 2) can all be expressed 
in terms of n and A;nb- To swiftly outline the general features of the observed charged multiplicity 
distributions it is appropriate to look at the most recent measurements and, in particular, at the 
results obtained at the LHC [23, 24, 25]. Consider in particular the data 4 of the Alice collaboration 
[23] which are summarized in Tab. 1 where the obtained results are also compared with the data of 
the UA5 experiment [26] for pp. In Tab. 1 (iV c h) coincides with n of Eq. (26) (only central values 
are reported). By looking at Tab. 1 three important points should be noticed. For fixed central 



4 We shall not dwell here on the details of the various experimental analyses such as the occurrence that the maximal 
multiplicity is always finite (i.e. n never goes to infinity) and, in this sense, all the charged multiplicity distributions 
reconstructed from the experimental data should belong to a truncated set of multiplicities for fixed energy and in central 
(pseudo)rapidity intervals. 



Data 


UA5 pp ^/s = 0.9 TeV 


Alice pp y/s = 0.9 TeV 


Alice pp y/s = 2.36 TeV 


(N ch ) for \r]\ < 0.5 
&NB for 1 77 1 < 0.5 


3.61 

1.50 


3.60 
1.46 


4.47 
1.25 


(iV ch ) for \r,\ < 1.0 

&NB for \l]\ < 1-0 


7.38 
1.62 


7.38 
1.57 


9.08 
1.37 


<iV ch ) for |r?| < 1.3 
/cnb for |t?| < 1-3 


no data 


9.73 
1.67 


11.86 
1.41 



Table 1: The parameters of the negative binomial distribution as they arise in different central 
(pseudo) rapidity intervals and different energies for the UA5 experiment and for the Alice experi- 
ment. The data are the ones published in [23] (see in particular the second paper). 



(pseudo)rapidity intervals, /cnb decreases as the centre of mass energy y/s increases while, as well 
known, (N&) increases (logarithmically) with y/s. For fixed centre of mass energy y/s, /cnb and (AT ch ) 
are both increasing with \n\, i.e. the absolute value of the (pseudo) rapidity 77. The two preceding results 
are consistent with the trends already observed in the GeV energy domain (see, e.g. [19, 20, 21]). 

When &nb decreases the tail of the distribution increases (i.e. correlations are present for large 
multiplicities in comparison with the Poisson case). There is a fourth point which is relevant for the 
forthcoming considerations: in the case of heavy ion collisions (see, e.g. [27] for the cases of Au+Au 
and Cu + Cu up to energies y/s = 200 GeV) /cnb is quite large (i.e. 1/^nb — > 0, and, according to Eq. 
(26), g<?) ->• 1). 

By comparing the multiplicity distributions in strong interactions with their gravitational coun- 
terpart various interesting connections emerge and will now be swiftly summarized. The inequality 
established in Eq. (25) for of gravitational multiplicity distributions holds also for hadronic multiplic- 
ity distributions in strong interactions (see Eq. (26)); in both cases the Bose-Einstein distribution (i.e. 



(2) 



2) and the Poisson distributions (i.e. g 



12) 



1) are, respectively, the signatures of maximal 



9 

and minimal correlation among the produced particles. The maximally correlated case in gravitational 
multiplicity distributions does not necessarily entail the presence of a mixed state but can even be 
the result of a pure state with zero charge. For hadronic multiplicity distributions, the decrease of 
/cnb (when the centre of mass energy of the collision increases and the central rapidity intervals are 
fixed) is pointing towards the Bose-Einstein limit (i.e. /cnb — > !)• For the gravitational multiplicity 
distributions the minimally correlated case corresponds to the situation where the number of species 
present in the initial state dominates over the average number of produced pairs. For hadronic mul- 
tiplicity distributions the minimally correlated case is expected in heavy ion collisions when /cnb is 
pretty large and the dispersion of the distribution roughly coincides with the average multiplicity 



(Poisson limit). While in the hadronic case &nb can a l so be non-integer, in the gravitational case /cnb 
assumes only integer values and counts the charges of the initial state: this is the main difference 
between the hadronic and the gravitational multiplicity distributions. It is nonetheless remarkable 
that the distorted negative binomial distribution (deduced in Eq. (15)) has a slight excess for low 
multiplicities while it is comparatively smaller than the standard negative binomial distribution for 
large multiplicities. The gravitational multiplicity distributions, as their hadronic counterparts, all fall 
within the class of infinitely divisible distributions and they can be, in turn, classified in terms of the 
matrix elements of the positive (dscrete) series of SU(1,1). In the case of gravitational multiplicity 
distributions there is an absolute asymptote for 1/&NB, i- e - 

J_^l, c 2 ^2 + j±-, (27) 

where C 2 is the second normalized moment (see Eq. (26) and discussion thereafter). Assuming that 
the trend established at the LHC will be confirmed, it would be interesting to understand if such an 
asymptote is also present for hadronic multiplicity distributions (in central (pseudo)rapidity intervals) 
as the centre of mass energy in pp collisions increases. 
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